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Abstract 

In a recent paper (Tran et al., Ann. Phys. 311(2004)204), some asymptotic number theoretical 
results on the partitioning of an integer were derived exploiting its connection to the quantum 
density of states of a many-particle system. We generalise these results to obtain an asymptotic 
formula for the restricted or coloured partitions p s k {n), which is the number of partitions of an 
integer n into the summand of s th powers of integers such that each power of a given integer may 
occur utmost k times. While the method is not rigorous, it reproduces the well known asymptotic 
results for s = 1 apart from yielding more general results for arbitrary values of s. 
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I. INTRODUCTION 



It is well known that in one- dimension, the density of states of a system of ideal bosons 
confined in a harmonic trap at a given energy E, in the N — > oo limit, is the same as the 
number of ways of partitioning an integer n = E into a sum of integers The result is the 
celebrated Hardy- Ramanuj an formula|2j. An in-depth analysis of the fundamental additive 
decomposition of positive integers by sums of other positive integers with various restrictions 
is given in Ref.(|3|). 

n 

In a recent paper Tran et al |4J exploited the connection between quantum many particle 
density of states and integer partitions in number theory to obtain more general results for 
the additive decomposition or partitioning of an integer. It was shown that the asymptotic 
density of states for ideal bosons in a system with a power law energy spectrum, is analogous 
to p s {n), that is, the number of ways of partitioning an integer n as a sum of s powers of 
integers. It was further shown that distinct partitions d s (n) correspond to the asymptotic 
density of states of a pseudo fermion-like system in a power law energy spectrum. The name 
pseudo fermion is coined to indicate that the system's density of states was constructed by 
applying the exclusion principle only to the particle distribution. The hole distribution was 
neglected while allowing any number of particles in the ground state, i.e., a collapse of the 
ground state, where the Fermi energy Ep is set to zero. For the rest of the paper we follow 
the same procedure as in Q] where only particle distribution, relevant to partitions, is taken 
into account but not the hole distribution. We note that the method is not rigorous but 
yields correct asymptotic results in the well known cases. 

In this paper we extend the analysis to obtain an asymptotic formula for restricted or 
coloured partitions. We define p s k {n) as the number of partitions of an integer n into a sum 
of s th powers of integers, where each power of a given number occurs utmost k times. The 
method used here is similar to the one used in the previous paper J^J. The leading exponential 
approximation for the asymptotic partition of an integer into integer summands, s — 1, has 

n 

been discussed recently by Blencowe and Koshnick[5j. 

We first construct the partition function for a statistical system in which the maximal 
occupancy of each state is given by k. Obviously k = 1 corresponds to a fermionic system 
where as there are no restrictions on k for a bosonic system. Indeed the partition function 
obtained by this minimal restriction is the same as that obtained by Gentile [6] and the 
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corresponding statistics is know as Gentile Statistics. The partition function of Gentile 
statistics also has the property that it nicely interpolates between the Fermi and Bose 
Statistics. While the asymptotic formula, obtained using Gentile statistics, reduces to the 
result for distinct partitions in the fermionic, k — 1, limit. The unrestricted or bosonic 
partitions are obtained by taking the k — » oo limit. 

II. MANY-PARTICLE DENSITY OF STATES 

Before we proceed with the main theme of this paper, we recall briefly results from the 
previous analysis for the sake of completeness. Consider a system of N particles, where N 
is very large. In general, the particles may obey any arbitrary statistics. The canonical 
iV— particle partition function is given by 



where (3 is the inverse temperature, E\ are the eigenenergies of the N— particle system and 
Pn(E) = ^2i$(E — e\ NS) ) is the N— particle density of states. The density of states Pn(E) 
may therefore be obtained through the inverse Laplace transform of the canonical partition 
function 



In general, it is not always possible to do this inversion analytically. As in the case of 
single-particle density of states, the many particle density of states may be decomposed 
into an average (smooth) part, and an oscillating part jjj. The smooth part ~p N (E) may 
be obtained by evaluating Eq.(J2J) using the saddle-point method pf. However, unlike the 
one-particle case, where the oscillating part may be obtained using the periodic orbits in a 
"trace formula" |?|, it remains a challenging task to find an expression for the oscillating 
part Spn{E) In what follows we shall use the saddle-point method to obtain the smooth 
part of the density of states, ~p N (E) and identify it with the number of ways the energy E 
is partitioned, on the average, among N particles. 
Defining the entropy, S(/3) as 




(1) 




(2) 



S{(3)= (3E + \ogZ N . 



(3) 
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and expanding the entropy around the stationary point (3q and retaining only up to the 
quadratic term in Eq.(J2J yields the standard result ^] 

= ^PiS(Po)] (4) 
where the prime denotes differentiation with respect to inverse temperature and 

The above analysis holds for any N particle system irrespective of the nature of the single 
particle spectrum. Here after we omit the the subscript N since we are only interested in 
the asymptotic limit iV — > oo. 



III. RESTRICTED PARTITIONS 



We consider a statistical system in which a given quantum state has a maximal occupancy 
of k in each single particle state. Thus when k — 1 the system corresponds to a system of 
fermions whereas there are no such restrictions on k for a system of bosons except for 
symmetry considerations. We now construct the partition function for such a system which 
also incorporates the property of interpolation between Bose and Fermi statistics as k takes 
on different values. 

In order to relate to the problem of integer partitions, we consider a system with a single 
particle spectrum given by e m = m s , where the integer m > 1, and s > 0. The energy 
is rendered dimensionless by appropriate choice of units. For example, when s = 1 the 
spectrum can be mapped on to the spectrum of a one dimensional oscillator, where the 
energy is measured in units of hio. For s — 2, it is equivalent to setting energy unit as 
h 2 /2M, where M is the particle mass in a one dimensional square well with unit length. 
These are the only two physically interesting cases. As in |4|, we however keep s arbitrary 
even though for s > 2 there are no quadratic Hamiltonian systems. 

We define a Restricted or Coloured Partition of an integer n as its additive decomposition 
into integer powers where no one integer is repeated more than k times. For example consider 
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partitions of an integer, say N = 6 as a summand of integers: 



6, 5 + 1, 4 + 2, 3 + 2 + 1, 


(Pi(6) 


= 4) 


3 + 3, 4 + 1 + 1, 2 + 2 + 1 + 1, 


04(6) 


= 7) 


3 + 1 + 1 + 1, 2 + 2 + 2, 


04(6) 


= 9) 


2 + 1 + 1 + 1 + 1, 


(ri(6) 


= 10) 


1 + 1 + 1 + 1 + 1 + 1, 


(pj(6) 


= H) 



(6) 

We are interested in the relation between the asymptotic expression for such restricted 
partitions and the density of states of the system which allows utmost k occupancy of the 
single particle states. 

The relevant physical system for obtaining partitions is thus one in which at zero temper- 
ature all the particles are in the ground state whose energy is set to zero. At any excitation 
energy, particles are excited from the ground state and distributed in the single particle 
states with maximal occupancy given by k in each state. The number of ways in which such 
a distribution can be achieved is the density of states at that energy (as also the number of 
ways of partitioning the energy). 

The partition function for such a system in the limit of the number of particles iV — ► oo 
is given by 

oo 

Zoo(P) = Y[[l + exp(-(3m s ) + --- + exp(-kPm s )} 

m=l 

oo k 

= n E ex p(-^ ms )' ( ? ) 

m=l n=0 

where the single particle spectrum is given by a power law. We note that the partition func- 
tion in Eq.(J2J) is indeed the grand-canonical partition function with the chemical potential 
fi = 0. In particular k = 1, oo correspond to the two cases considered in the earlier paper j^j. 
Note that with k arbitrary this is indeed the partition function for the well-studied Gentile 
statistics 0. 

By setting x = exp(— (3), the partition function may be written as 

ZM = J2pl(n)x" = H [ J . (8) 

n=l n=l L J 

In the asymptotic limit of large number of particles, p|(n) is the number of ways of parti- 
tioning n. 



Using Eq.QjWe obtain 

oo oo 

S = (3E + J2 ln[l-exp(-(A; + l)/3n s )]-^ ln[l - exp{-(3n s )]. (9) 

n=l n=l 

We now evaluate the sums approximately using the Euler-Maclaurin series. We have 

oo oo 

- H 1 - exp(-/3n s )] = -J2 H 1 ~ exp(-/3(n + l) s )] 

n=l n=0 

/oo i °° B 

ln[l - exp(-/?(x + l) s )}dx - - ln[l - exp(-/3)] + s £ ^^y, (10) 

where i?2fe are the Bernoulli numbers. The integral in the Euler-Maclaurin expansion may 
be evaluated as follows: 

POO 

- / ln[l - exp(-/5(x + l) s )}dx 
Jo 

roc pi 

— — ln[l - exp(-[3x s )}dx + / ln[l - exp(-(3x s )]dx 

Jo Jo 



where 



c( a ) = r(i + i)c(i + i) (12) 



given in terms of the Riemann zeta function. In the high temperature limit, that is (3 — > 0, 
we have 

~y ln[l - exp(-f3(x + l) s )]dx « ^ + ln[l - exp(-/3)] - s + C(/3) (13) 

The series involving Bernoulli numbers may be evaluated approximately as follows: The 
Euler-Maclauring series applied to logarithm of Gamma function gives, 

1 00 B 

ln(T(n + 1)) = n ln(n) - n + - ln(27r) + £ 2fc(2fc Z^j^n 

We may consider this as the expansion of zero by putting n — 1 and get 

y 7 B ? k - = s--ln(27r). (15) 



(14) 



00 

s 

k 



Thus keeping terms up to 0(/3), and using the above arguments the last term in Eq.© 
may be written as 

- f> - expH*,-)] * + Nl-^-g)] _ I ln(2x) + om (16) 

n=l 



One may also get the above result directly using symbolic math program MAPLE. The infi- 
nite sum in the second term of Eq. @ is also evaluated similarly by scaling /3 appropriately. 
Note that the constant term does not get scaled. 

Thus the entropy to the leading order in (3 is given by 

S = (3E + - i ln[l - <Kp(-(* + l)(3)} + ~ ln[l - ex P (-/3)] + 0((3) , (17) 

where 

a = l — r . (18) 

(Jfe + l)* 

We note that the expansion in powers of (3 is equivalent to taking the high-temperature 
limit though this limit should be employed cautiously when j3 is weighted by k which can 
in principle take large values. The last term of 0(f3) in the RHS of Eq. (jl7j) causes a shift 
in the energy E by a constant (for example 1/24 in the s=l case) in S. In the asymptotic 
formula for the level density for large E, obtained by the saddle point method, this may 
be ignored. Note, however, that such a shift was included by Rosenzweig[10] to improve 
the Bethe formula jlll] for the nuclear level density. We further note that the k — *> oo, or 
the bosonic limit, is special since this limit involves precisely the sum given in Eq. (jl6|) and 
not the difference of two series as in Eq. As a result the (3 independent term, namely 
— (s/2) ln(27r), appears only in this limit and not for any finite k. Indeed, as we shall 
see below this term is crucial in getting the prefactor correctly in the Hardy-Ramanujan 
Formula. 

The saddle point on the real axis [3] is evaluated by taking the derivative of the entropy, 

s '^ = E -\w&' (19) 

where only the leading terms in (3 are retained. Equating the derivative to zero, the saddle- 
point in (3 is given by 

A-fsW** =«.*-*, (20) 



where 



sE 



n s = ) . (21) 



s 

Substituting this in the saddle point expression for the density of states in Eq.(jU, the 
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asymptotic density of states is given by 



exp 



K s (s + l)Ei+t 



22) 



y/2n(s + i)£(3s+i)/(*+i)[i - exp(-(fc + 1)k s E- s /(«+i))] 
where the subscript k in p indicates the main property of the statistical system under consid- 
eration. We have kept the exponential term under the square-root sign in the denominator 
to indicate the interpolation property of the density of states between the Bose (k — ► oo) 
and the Fermi (k = 1) limits even though keeping the exponential may not be consistent 
with the order of the expansion. The correct expression is given below. 

For finite k there always exists an energy E which is large enough such that the following 
approximation is useful: 

kJs + 1)E^ 



exp 



(23) 



v /2vr(s + l)(A; + l)E( 2s + 1 )/( s + 1 ) ' 
Note that the above expression is .den.ee! t o the E,(23) for «.(„) in Q eo rre spo ndillg to 
the special case of k = 1. The degeneracy or k dependence in the above equation is also 
hidden in the parameter k s . An asymptotic formula derived in Ref.(^]) based on Gentile 
statistics for the special case of s = 1 is the same as the one given in Eq.(|23|). However, to 
the best of our knowledge, no general formula for arbitrary s and k exists. 

In particular for s = 1, k = 1 the above equation reduces to the well known formula P.ll2^ 
for distinct partitions of an integer into a set of integers 

exp 7TW | 

p{E) « d(n = E) = 4 - 3 L 1/4 y 4 • (24) 

The bosonic limit may be obtained from Eq.()22j) by taking the limit k — > oo. However 
we also note that this limit should be taken in such a way that the (3 independent factor in 
Eq.(|16|) is also included. With this proviso we obtain 

k s (s + l)E^i 



pI s \e) 



exp 



(25) 



which is the celebrated Hardy- Ramanujan asymptotic formula for integer partitions. In 
particular for s = 1, we have the well known result (see for example Ref.(|l2|). 



p(E) w p(n = E) 



exp 



7T 



aVze 



(26) 



In Fig. (^) we show a comparison between the exact pk{n) (continuous line), and the 
asymptotic density of states, p k (n) (dashed line) obtained from Eq. (J23|) for various values 
of k when the power s — 1. For smaller k the asymptotics is reached earlier than for larger k 
as can be seen easily from the figures. In fact for k = 1 the formula is almost exact even for 
n small. We note that the numerical estimates for the density based on Eq. (j22j) are usually 
higher than the exact values and asymptotically reach the exact value later than the density 
calculated based on Eq. (j23j) which is shown in the figure. The saddle point evaluated by 
keeping only the two leading terms in (3 ensures that the Eq. ()23|) is more accurate for small 
values of k though asymptotically both equations yield the same result. 

Similarly, in Fig. (j2J), the computed p\{n) is compared with ~p\{n). We note that in this 
case the fluctuation about the mean is larger for small values of k. In general the fluctuations 
decrease with increasing k for a given s but increase with increasing s for a given k as shown 
in FigJ3] where we have plotted the asymptotic density of states for s = 3 and various values 
of k. 

IV. SUMMARY 

To summarise, this paper extends the results of Ref. Q] to the case of restricted or 
coloured partitions. We identify the relevant generating function for coloured partitions as 
the partition function of gentile statistics. The resulting asymptotic or smooth density states 
interpolates between the Bose (arbitrary partitions) and Fermi (distinct partitions) systems 
as a function of the maximal occupancy of the quantum state. 
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thank the hospitality of the Institute of Mathematical Sciences where this work was done. 
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FIG. 1: Comparison of the exact p\{n) (solid line) and the asymptotic p\{n) (dashed line) for s 
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